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Abstract 

The pedagogic two stste system of the ammonia molecule is used to illustrate the phenomenon 
of environment induced molecular localization in pyramidal molecules. A semiclassical model is 
used to describe a gas of pyramidal molecules interacting via hard ball collisions. This modifies 
the tunnelling dynamics between the classical equilibrium configurations of an isolated molecule. 
For sufficiently high pressures, the model explains molecular localization in these classical configu- 
rations. The decrease in the inversion line frequency of ammonia, noted upon increase in pressure, 
is also studied. 

PACS numbers: 03.65.Xp, 33.20.Bx, 73.43.Nq, 33.55.Ad 
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I. INTRODUCTION 



An isolated ammonia molecule would periodically switch between its classical equilibrium 
states [lj]. Its stationary Eigenstates are delocalized symmetric and antisymmetric ones 0. 
The energy difference between these states account for the absorption frequency of NH 3 at 
low pressure. It is the ideal system to teach the quantum mechanics of two state systems. 

But, experimentally it has been observed that, as pressure is raised, the inversion line 
shifts and ultimately dissappears. The molecule localizes in one of the classicaly stable 
configurations. This phenomenon has been investigated for a long time 0] . 

A qualitative description using dipole-dipole interaction had been suggested Q. It can 
be quantitatively explained using the Keesom Energy 0], if the interaction terms in a 
mean field Hamiltonian cause the molecule to be stationary in one of the classical localized 
configurations. A lot of work [3] [1] [s| has has gone into explaining the phenomenon. The 
role of collisions in inducing molecular localization has also been suggested [hJ. Perturbative 
terms have also been used to represent weak collisions [Tlj ] . 

The purpose of this work is make the students, with basic knowledge of quantum me- 
chanics, familiar with this active branch of research. A semiclassical model is presented for 
the molecules interacting via simple collisions in a Van der Waals' Gas, hence introducing a 
mean field interaction constant, representing hard ball collisions. The resulting Hamiltonian 
is used to find out the new eigenstates as well as the inversion line frequency. 

II. ISOLATED MOLECULE 

A. Effective Potential 

For an isolated molecule the one dimensional motion of the nucleus N across the plane 
containing the three H nuclei can be considered separately from the other degrees of freedom 

Q. 

The effective potential energy 0] is a function of the distance x between the N nucleus 
and the plane containing the H nuclei. The symmetry of the system about x = dictates 
V(x) to be an even function of x. V(x) has two minima corresponding to its classically 
stable configurations. Let us assign this minimum value as zero. 

The N atom is repelled by the H atoms if it approaches their plane. This is represented 
by an energy maxima of height V\ at x — 0. The chemical bonding force ensures an increase 
in V beyond the two minima. Hence, forming a double well potential. 

B. Hilbert Space 

We say that the system is in state \L) when it is on one side and \R) when it is on the 
other side. 




(1) 



(2) 



Hence, forming the basis vectors of a two dimensional Hilbert space. 
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C. Hamiltonian 



We choose the Hamiltonian operator on the two-state Hilbert space to be 



AE 



.X 



(3) 
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where a x is the Pauli matrix 




1 

1 



) 
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D. Eigenstates 



The Eigenvectors corresponding to the Hamiltonian formulated in Eq. ([3]) are 




(6) 



(5) 



with Eigenvalues —AE/2 and +AE/2, respectively. 
E. Absorption Frequency 

The inversion spectra of Ammonia is thus explained in terms of the energy difference 
between these two stationary states. 



III. HARD BALL COLLISIONS 

The gas is modelled as a system of nearly independent molecules, such that each molecule 
is subjected to an external field representing random hard ball collisions with the rest of the 
molecules. This interaction is to be determined self-consistently. Then the linear response of 
this model to an electromagnetic perturbation is to be studied to obtain the absorption fre- 
quency. The inter-molecular interaction is found to control this frequency, hence explaining 
the shift in the the inversion line. 

The random motion of the molecules are faster than the inversion times, so on the time 
scales of the inversion dynamics the molecules feel an effective interaction arising out of the 
random collisions with the other molecules. 

In the representation chosen for the Pauli matrices, the localizing effect of the collision 
between two molecules % and j can be represented by an interaction term of the form ofcrj, 
where a z has localized Eigenstates 




(9) 



(8) 
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IV. INTERACTION ENERGY 



For extremely rarefied gases, the molecules are considered to be negligibly small in size 
and the gas is supposed to have only kinetic energy and no potential energy. The equation 
of state for 1 mole of such a gas is 

PV = RT, (10) 

where P is the pressure and V is the molar volume. But, the finite size of the gas molecules 
effectively reduce the volume available for the free motion of the molecules. If the molecules 
have a diameter d, no other molecular center of mass can come within a sphere of radius 
d, centered at the center of mass of the molecule under consideration. According to Van 
der Waals, considering the molecules to be classical mechanical hard balls which only app ly 
contact forces and do not interact when away from each other, the equation of state [l2j of 
the gas hence becomes 

P(V-b) = RT, (11) 
where b is defined as half the excluded volume, due to the finite size of the molecules. Hence, 

LA 

b=--nd 1 d 2 d 3 , (12) 

where the di are the principle collision diameters and L, the Avogadro's number, is the 
number of molecules present in our system. 
Expanding Eq. ffTTl) we have 

PV-Pb = RT, (13) 

Here, Pb is the energy excluded from the equation of state due to hard ball collisions, so, it 
must be accounted for. If C be the mean excluded pressure energy per molecule, then 

LC = Pb. (14) 

Substituting, b from Eq. (fT2l . we have 

PA 

C = —-7idid 2 d 3 . (15) 

Since, we are considering only one degree of freedom for the molecule, d\ = d 2 = d, while 
the other principle collision diameter d 3 is to be determined self consistently from the state 
of the given molecule. Now, the molecule has a maximum size for the localized Eigenstates 
of a z . So we write, 

d 3 = (i/;\<7 z \ rp)d (16) 

Hence, the mean field interaction energy, per molecule, arising out of the hard ball collisions, 
becomes 

C = C Q>\o*\il>) i (17) 
where, the mean field interaction constant 

Co = (18) 
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V. MEAN FIELD HAMILTON! AN 



In the mean-field approximation laid out in Eq ffTTj) . we obtain the total Hamiltonian 

AE 

h(ip) = — a x - Ca z . (19) 

Substituting C from Eq. (fT8l) . we have 

hW = -^-o» ~ C (4j I a* | ^)a\ (20) 

where is the state of the molecule being considered and C is given by Eq. (TTS1) . The state 
is to be determined self-consistently by solving the nonlinear Eigenvalue problem associated 
with this Hamiltonian. 



VI. EIGENSTATES 

The solution of the Eigenvalue problem associated with the Hamiltonian gives the fol- 
lowing results. 

If Co < AE/2, then 

|Ai> = |1>, (21) 
the symmetric Eigenfunction, with Eigenvalue /ii = — 4p, is the only ground state. While, 

|A 2 ) = |2>, (22) 

with Eigenvalue /i 2 = 4p forms the first excited state. 

But, if C > AE/2, then there are two degenerate ground states, 



,1 AE . . 1 AE . . 
|A^> = A /- + ^|l> + ^- - _ |2> (23) 



and 



with /i 3 = /i 4 = — C . 

Keeping this dependence on Co in mind, the critical value 

C cr = AE/2 (25) 

is defined. This distinguishes the following two cases. For Co < C™, the ground state of the 
system is the delocalized symmetric state of an isolated molecule. For Co > C^, there are 
two different states which form the stationary states of the system. 

These states described in Eqs. ( 123H2^|) . transform into one another under the action of 
the inversion operator a x , 

|Af) = ^|A 4 L ) (26) 

hence called chiral. 
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When Cq ^> Cq, these new ground states become localized, 

lim |A£) = \L) lim I A?) = \R) (27) 

AE/C -*Q AE/Co^O 1 4/ 1 ' v ' 

The emergence of the states described in Eqs. (I23H24P imply a bifurcation of the ground 
state at a critical interaction Co = Cq . Using the expression for Co from Eq. (TT8]) . this 
bifurcation can be obtained above the critical pressure 

Using d = .274 nm, AE corresponding to 23.8 GHz, we have P cr = 1.8 atm. 



VII. INVERSION LINE FREQUENCY 

Electro-magnetic radiation of angular frequency ujq adds a perturbation term j^] to the 
Hamiltonian described in Eq. (1201) 

h em (t) = ef(t)a z (29) 

where e is a small parameter and f(t) = 9(t) cos(u t), 6{t) being the Heaviside function. 
The choice of a dipole coupling approximation, h em oc a z , is justified for a radiation of 
wavelength long with respect to the molecular size. 

Under the effect of the perturbation the isolated molecule state |A(i)) evolves according 
to the time-dependent nonlinear Schrodinger equation 

it^^ = [ h (X(t))+ef(t)a^\X(t)) (30) 

where h(X) is the Hamiltonian for a molecule in the mean field described in Eq. ( 1201) . 

The linear response to the perturbation is expressed by the generalized susceptibility 
1Z{uj) = Si(u) / f (u) , where f(w) and S\(u) are the Fourier transforms of f(t) and <5>i(i), 
with Si(t) defined by 

S(t) = (X(t) I a* | A(t)> 

= S (t)+eS 1 (t) + ... (31) 

Let us assume that at time t — each molecule is in the delocalized ground state |Aq) = 
The solution of the dynamical Eq. (|30l) . with the initial condition |A(0)) = |1), gives 

. 2AE . 

U{UJ) = (M 2 -(A^-2C AE)- (32) 

The generalized susceptibility has a unique pole at positive frequency which corresponds to 
the inversion line frequency 



h V AE , 

Substituting the value of Cq from Eq. (ITBl) this equation becomes 



which compares well with experimental data [13| [14 
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VIII. CONCLUSIONS 



Students are customarily exposed to the symmetric and antisymmetric ground states 
of a system. Quantum mechanics is taught as the basis of molecular structure. Yet a 
large class of molecules are found to be chiral, not symmetric. It is beyond the scope of 
undergraduate syllabii, as well as this work, to address the quantum origins of chirality 
of organic molecules. But, it can be instructive for students to learn about environment 
induced localizations through the use of simple models such as this. 
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